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UNCLASSIFIED  ABSTRACT 


BUCKLING  OF  CIRCULAR  CYLINDRICAL  TR-0158(S3820-10}-1 

SHELLS  WITH  MULTIPLE  ORTHOTROPIC  September  1967 

LAYERS  AND  ECCENTRIC  STIFFENERS, 
by  Robert  M.  Jones 

An  exact  solution  is  derived  for  the  buckling  of  a  circular  cylindrical 
shell  with  multiple  orthotropic  layers  and  eccentric  stiffeners  under 
axial  compression,  lateral  pressure,  or  any  combination  thereof. 
Classical  stability  Mieory  (membrane  prebuckled  shape)  is  (..sed  for 
simply  supported  edgs  boundary  conditions.  The  present  theory 
enables  the  study  of  coupling  between  bending  and  extension  due  to 
the  presence  of  different  layers  in  the  shell  and  to  the  presence  of 
eccentric  stiffeners.  Previous  approaches  to  stiffened  multilayered 
shells  are  shown  to  be  erratic  in  the  prediction  of  buckling  results 
due  to  neglect  of  coupling  between  bending  and  extension. 
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NOMENCLATURE 

=  ring  spacing  (Figure  1) 

=  cross-sectional  area  of  a  stiffener 
=  coefficients  in  stability  criterion  [Eq.  (18)] 

=  stringer  spacing  (Figure  1) 

=  extensional  stiffness  of  the  layered  shell 
=  extensional  stiffness  of  the  orthotropic  stiffness 
layer  in  the  x-(y-)  direction 

=  in-plane  shearing  stiffness  of  the  orthotropic  stiffness 
layer 

=  coupling  stiffness  of  the  layered  shell 
=  bending  stiffness  of  the  layered  shell 
=  bending  stiffness  of  the  orthotropic  stiffness  layer 
in  the  x-(y-)  direction 

=  twisting  stiffness  of  the  orthotropic  stiffness  layer 

=  Young's  modulus  of  a  stiffener 

=  Young's  moduli  in  x  and  y  directions, 

til 

respectively,  of  the  k  ‘  shell  layer 

=  shearing  modulus,  E/(2(l  +  v)},  of  a  stiffener 

th 

=  shearing  modulus  of  the  k  shell  layer  in  x-y  plane 
=  moment  of  inertia  ot  a  stiffener  about  its  centroid 
=  torsional  constant  of  a  stiffener 


A  comma  indicates  partial  differentiation  with  respect  to  the  subscript 
following  the  comma.  The  prefix  6  denotes  the  variation  during 
buckling  of  the  symbol  which  follows. 
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NOMENCLATURE  (Continued) 


j<k  =  function  of  material  properties  of  the  k 

ij 

layer  [Eq.  (2)] 

L  =  length  of  circular  cylindrical  shell  (Figure  1) 

m  =  number  of  axial  buckle  halfwaves 


M  ,  M  , 
x’  y’ 


M  ,  M 
xy’  yx 


N  ,  N  ,  N 
x’  y  xy 

N  ,  N 
x’  y 


u,  v,  w 


x,  y,  -a 


(x’  cy’  ^ xy 

V  V  (3 


=  moments  per  unit  length 

=  number  of  circumferential  buckle  waves 
=  number  of  layers 
=  in-plane  forces  per  unit  length 

=  applied  axial  and  circumferential  forces  per  unit 
length 

=  external  or  hydrostatic  pressure 
=  shell  reference  surface  radius  (Figures  1  and  2) 

=  thickness  of  k^  shell  layer 

-  axial,  circumferential,  and  radial  displacements  from 
a  membrane  prebuckled  shape 

=  axial,  circumferential,  and  radial  coordinates  on  shell 
reference  surface  (Figure  I) 

-  distance  from  stiffener  centroid  to  shell  reference 
surface  (Figure  1),  positive  when  stiffener  on  outside 

-  strains 

-  variations  in  reference  surface  strains  [Eq.  (5)] 

distance  from  inner  surface  of  layered  shell  to  outer 
th 

su  rfac.e  of  k  la  ye  r 

distance  from  inner  surface  of  layered  shell  to 
reference  surface 
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NOMENCLATURE  {Continued) 


) 


VxyB*VyxB) 


=  Poisson's  ratio  for  contraction  in  the  y(x)  direction 
due  to  tension  in  the  x(y)  direction 
=  so-called  extensional  Poisson's  ratio  for  contraction 
in  the  y-(x-)  direction  due  to  tension  in  the  x-(y-) 
direction 

=  so-called  bending  Poisson's  ratio  for  curvature  in  the 
y-(x-)  direction  due  to  moment  in  the  x-(y-)  direction 


<r  ,  a-  ,  t  =  stresses 

x’  y’  xy 

^1’  ^2’  *3  =  var*ati°ns  reference  surface  curvatures  [Eq.  (6)] 


Superscript 

k 

=  k1*1  shell  layer 

Subscripts 

k 

=  k^  shell  layer 

r 

=  ring 

s 

=  stringer 
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SECTION  I 


INTRODUCTION 

The  first  work  in  the  area  of  stability  of  eccentrically  stiffened 
shells  was  done  by  Van  der  Neut  (Ref.  1)  about  twenty  years  ago. 
However,  his  conclusion  that  the  buckling  load  under  axial  compres¬ 
sion  of  an  externally  stiffened  shell  can  be  as  high  as  two  or  three 
times  that  of  an  internally  stiffened  shell  went  essentially  unnoticed. 
More  recently,  Baruch  and  Singer  (Ref.  2)  and  Block,  Card,  and 
Mikulas  (Ref.  3)  presented  theories  which  are  considered  basic  in  the 
field.  Since  1965,  work  in  the  area  of  eccentrically  stiffened  shells 
has  expanded  so  much  that  it  is  impractical  to  mention  more  than  a 
few  significant  papers.  McElman,  Mikulas,  and  Stein  (Ref.  4) 
extended  the  original  work  to  include  the  effect  of  stiffeners  on  vibra¬ 
tion  and  flutter.  Correlation  between  theory  and  experiment  was 
reported  for  static  buckling  loads  by  Card  and  Jones  (Ref.  5).  The 
effect  of  initial  imperfections  was  considered  by  Hutchinson  and 
Amazigo  (Ref.  6).  Block  (Ref.  7)  treated  discrete  ring  spacing, 
prebuckling  deformation,  and  load  eccentricity.  Finally,  plastic 
buckling  was  discussed  by  Jones  (Ref.  8). 

The  object  of  the  present  paper  is  to  extend  previous  theories 
to  consideration  of  stability  of  circular  cylindrical  shells  with 
multiple  orthotropic  layers  and  eccentric  stiffeners  (see  Figure  1). 
Classical  stability  theory,  which  implies  a  membrane  prebuckled 
shape,  is  used  for  the  simply  supported  edge  boundary  conditions 
6Nx  v  -  w  =  6Mx  =  0.  The  layers  have  orthotropic  material 


1 


properties  with  the  principal  axes  of  orthotropy  coincident  with  the 
shell  coordinate  directions.  In  accordance  with  most  previous 
theories,  the  stiffeners  are  treated  as  isotropic  one-dimensional 
beam  elements  and  are  averaged  or  "smeared  out"  over  the  stiffener 
spacing.  The  torsional  rigidity  of  the  stiffeners  is  accounted  for  in 
an  approximate  manner.  The  present  theory  enables  the  study  of 
coupling  between  bending  and  extension  due  to  the  presence  of 
different  layers  in  the  shell  and  to  the  presence  of  eccentric 
stiffeners. 
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SECTION  II 


DERIVATION  OF  THEORY 

Expressions  are  obtained  for  the  variations  of  stresses  during 
buckling  in  the  k**1  layer  of  a  multilayered  shell  in  terms  of  the  variations 
of  strains  during  buckling.  Subsequently,  the  variations  of  stresses  are 
integrated  over  the  shell  and  stiffeners  in  order  to  obtain  expressions 
for  the  variations  in  forces  and  moments  during  buckling.  Finally,  the 
variations  in  forces  and  moments  are  substituted  in  Donnell-type 
stability  differential  equations  which  are  then  solved  to  yield  a  closed- 
form  stability  criterior  in  terms  of  the  geometric  and  material  properties 
of  the  stiffened  multilayered  circular  cylindrical  shell. 

A.  ORTHOTROPIC  STRESS -STRAIN  RELATIONS 

The  stress -strain  relations  for  an  orthotropic  material  can  be 
written  as 


<r  k  =  K  k  «  +  K.k_  « 

x  11  x  12  y 


„  k  T .  k  ,  „  k 
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y  12  x  22  y 
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xy  33  xy 
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wherein  the  superscript  k  denotes  the  k  layer.  The  quantity 


E  k  (E  k 

xx  yy 


is  Young's  modulus  in  the  x  (y)  direction,  Gk  is  the 

xy 
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»ic  !c 

plane,  and  v  ( v  )  is  the  Poisson's 
’  xy  yx 

ratio  for  contraction  in  the  y  (x)  direction  due  to  tension  in  the 

x  (y)  direction.  There  are  apparently  five  material  constants  per 

k  k 

layer;  however,  because  of  the  reciprocal  relations  (v  E^  = 
k  k 

Vyx  ®yy)»  there  are  actually  only  four  independent  constants. 

B.  VARIATIONS  OF  STRESSES  AND  STRAINS  DURING  BUCKLING 
During  buckling,  the  stresses  vary  from  their  prebuckling  values. 
Let  the  variation  be  denoted  by  6  ;  then,  from  Eq.  (1) 
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in  the  strains  during  buckling.  Because  of  the  Kirchhoff-Love  hypothesis, 
the  variations  in  strains  during  buckling  are 
6c  =  t .  +  z  X 
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6,y  =  <2  *  sX; 
6Yxy  =  <3  +  >*j 


(4) 


The  z  coordinate  is  measured  from  an  arbitrary  reference  surface 
(see  Figure  1 ).  In  Eq.  (4),  c^,  and  c^  are  the  variations  of  the 


reference  surface  strains 


u. 


<2  =  V’ 


u. 


+  w/R 


+  v, 


X 


(5) 
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and  X  j,  X^,  and  X^  are  the  variations  of  the  reference  surface 
curvatures 


X.  =  -  w, 

1  ’  xx 


X2  ’ 


w 


’  yy 


(6) 


X,  =  -  2w, 

3  ’  xy 


th 


Upon  substitution  of  Eq.  (4),  the  variations  in  stresses  in  the  k  layer 
can  be  written  as 


6<rxk  =  KU<'l+*Xl>  +  Kik2<‘2  +  zX2l 


6<rvk  =  Klk2<‘l+2Xl)  +  K2k2(<2  +  lX2> 


6rk  =  (£-  +zX,) 

xy  33  3  3 


(7) 


C.  VARIATIONS  OF  FORCES  AND  MOMENTS  DURING  BUCKLING 


The  variations  of  forces  and  moments  during  buckling  are  obtained 
by  integration  of  the  variations  of  stresses  over  the  shell  layers  and 
stiffeners.  The  effect  of  the  stiffeners  on  the  variations  of  forces  and 
moments  is  averaged  or  "smeared  out"  over  the  stiffener  spacing. 
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th 

where  denotes  the  thickness  of  the  k  layer  and  N  is  the  number 

of  layers.  The  variations  of  stresses  for  the  stiffeners  are  based  on 
uniaxial  isotropic  reductions  of  the  orthotropic  stress-strain  relations. 
The  integrations  in  Eqs.  (8)  and  (9)  yield 
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The  stiffnesses  in  Eq.  { 12)  are  due  to  Ambartsumyan  (Ref.  9)  and  depend 
on  the  location  of  the  reference  surface  (see  Figure  2).  The  reference 
surface  can  be  changed  by  varying  A  in  order  to  study  different 
loading  and  boundary  conditions.  Geier  (Ref.  10)  obtains  expressions  which 
arc  more  simple  in  appearance  than  Eq.  (1C),  but  which  arc  more 

difficult  to  utilize. 

LAYER 


Figure  2.  Cross  Section  of  an  N-Layered  Shell 
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D.  STABILITY  DIFFERENTIAL  EQUATIONS 

The  Donnell-type  stability  differential  equations  for  circular 
cylindrical  shells  subjected  to  combinations  of  axial  compression  and 
lateral  pressure  are 

(13) 

and  the  alternative  force  and  geometric  boundary  conditions  at  x  =  0 
and  L  are  chosen  from  the  following  sixteen  possibilities  (any  set  of 
four  alternatives  in  the  following  pairs  constitutes  a  set  of  boundary 


Upon  substitution  of  the  expressions  for  the  variations  of  forces  and 
moments  during  buckling  [Eqs.  (10)  and  (11)]  and  the  variations  of 
reference  surface  strains  and  curvatures  [Eqs.  (5)  and  (6)]  ,  the 
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stability  differential  equations  become 
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E.  STABILITY  CRITERION 


It  is  desired  to  find  the  solution  to  the  stability  differential 

equations  for  the  simply  supported  edge  boundary  conditions 

oN  -  /  -  w  -  6  M  -  0 

x  x 


ii 

£? 


1  i 


•  # 


The  following  buckling  displacements  satisfy  the  boundary  conditions  of 

Eq.  (16): 

u  =  u  cos{mrrx/L)  cos  (ny/R) 

v  =  v  sin(mTrx/L)  sin  (ny/R)  f  (17) 

w  =  w  sin(rmrx/L)  cos  (ny/R) 

(where  u,  v,  and  "w  are  the  amplitudes  of  the  buckling  displacements)  and 
are  substituted  in  the  stability  differential  equations  [Eq.  (15)3  •  1° 
order  to  obtain  a  nontrivial  solution  to  the  resulting  equations,  the 
determinant  of  the  coefficients  of  u,  v,  and  w  must  be  zero,  and  the 
following  stability  criterion  results: 
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=  (B1?/R)(rmr/L)  +  (Cn  +  ZSE  A  /b)(mir/L)‘ 

,2 


+  (Cu  +  2C33)(mTr/L)(n/R)‘ 


A,,  -  B._( rmr/L)2  +  (B„  +  E  A  /a)(n/R)2 
22  33  Li.  x  x 


2  1 


(C12  t  2C33)(rr.ir/Lr(n/R)  +  (1/R)(B,,  +  E_A_/a)(n/R) 


r  r 


(C,,  +  z  E  A  /a)(n/R): 
22  r  r  r 
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(18) 


(19) 

Cont. 

next 

paf,e 


A33  *  (DU  +  ESVb  *  S2EsAs/b)(mWL)4 

+  {4D33  +  2D12  +  GgJs/b  +GrJr/a)(m7r/L)2(n/R)2 

+  (D22  +  E  l  /a  +  zr2ErAr/a)(n/P.,4  +  (2C12/R)(rmr/L)2 

4-  (2/R)(C22  +  zrErAr/a)(n/R)2  +  (1/R2)(B22  +  E^/a) 

The  solution  represented  by  Eq.  (18)  reduces  to  the  solution  of 
Ref.  3  for  stiffened  single-layered  isotropic  circular  cylindrical  shells. 
In  addition,  stiffener  eccentricity  is  more  obviously  accounted  for  in  the 
foregoing  derivation  than  in  the  work  of  Geier  (Ref.  10) 

The  buckling  load  under  axial  compression  is  obtained  from 
Eq.  (18)  by  equating  N  to  zero  and  solving  for  N^.  Similarly,  the 
buckling  load  under  lateral  pressure  is  obtained  by  equating  to  zero 

and  solving  for  N^(N^  =  pR/t).  Finally,  the  buckling  load  under 
hydrostatic  pressure  is  obtained  by  equating  to  N^/2  and  solving  for 

N  .  In  addition,  if  N  (N  )  is  fixed,  the  critical  value  of  N  (N  )  can 
be  found.  In  this  manner,  an  interaction  curve  between  axial  compression 
and  lateral  pressure  can  be  obtained. 

Because  of  the  numerous  parameters  in  Eq.  (18)  and  the  need  to 
investigate  a  large  range  of  buckling  modes  to  determine  the  lowest 
buckling  load,  it  is  necessary  from  a  practical  standpoint  to  use  a  digital 
computer  for  numerical  work.  In  the  computer  program  (see  Appendixes 
A  and  C),  for  a  given  number  of  axial  halfwaves,  m,  and  circumferential 
waves,  n.  in  the  buckled  shape,  the  appropriate  buckline  load  is  found. 
The  number  n  is  varied  in  an  inner  DO  loop  for  a  fixed  m  until  ail 
relative  minima  of  the  buckling  load  are  found  within  a  given  rangi  cJ 


(19) 

Cont. 

from 

prev. 

page 


1  1 


values  of  n.  The  number  m  is  then  varied  in  an  cuter  DO  loop  so  that 
all  relative  minima  are  found.  Finally,  the  absolute  minimum  buckling 
load  is  selected  from  the  relative  minima. 
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SECTION  m 


NUMERICAL  EXAMPLE 

Because  of  the  many  geometrical  properties  in  the  theory, 
meaningful  general  results  cannot  be  presented.  Accordingly,  a  specific 
numerical  example  is  given  to  illustrate  application  of  the  theory.  The 
results  are  compared  with  results  of  previous  approaches  to  the  same 
problem. 

For  this  example,  the  stability  of  a  ring -stiffened  circular  cylindri¬ 
cal  shell  with  two  isotropic  layers  under  hydrostatic  pressure  is  considered. 
The  properties  of  the  layers  are 

Ej  =  44  x  106  psi  E2  =  2  x  106  psi 

vl  =  ^  V2  =  °.  4 

tj  =  0.  04  in.  t2  =  0.  3  in. 

The  rings  are  of  rectangular  cross  section  with  a  height  of  0. 25  inch 

and  a  thickness  of  0.  06  inch.  The  rings  are  on  the  inner  surface  of 
layer  one  and  have  the  same  material  properties  as  layer  one.  The  shell 
has  a  length  of  12  inches  and  a  radius  of  6  inches  to  the  middle  surface 
of  layer  one  (which,  in  this  case,  is  also  the  reference  surface). 

The  hydrostatic  buckling  pressure  of  the  above  configuration  is 
shown  as  the  solid  line  in  Figure  3  as  a  function  of  ring  spacing.  The 
results  shown  are  for  general  instability  (buckling  in  which  the  rings 
participate).  The  buckling  pressures  for  panel  instability  (buckling 
between  rings)  are  much  higher  than  the  present  results  and,  hence,  do 
not  govern  the  stability  of  the  present  configuration.  Other  failure  criteria, 
e.g.,  yielding,  are  ignored  for  the  purposes  of  this  illustration  of  the 
present  analysis  technique.  The  dashed  curve  in  Figure  3  represents 
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an  orthotropic  stiffness  approach  to  the  problem  and  is  from  3  to  9 
percent  lower  than  the  results  from  the  present  theory.  These  lower 
results  are  due  to  neglect  of  coupling  between  bending  and  exten¬ 
sion  of  the  layered  shell  and  the  eccentric  stiffeners  in  the  ortho¬ 
tropic  stiffness  approach.  The  solid  curve  with  a  single  dot 
represents  a  stiffened  shell  with  a  single  equivalent  Poisson's  ratio 
for  bending  ( =  0.  331)  used  in  both  layers  (Ref.  11)  and  is  from 
7  to  1 1  percent  lower  than  the  results  of  the  present  theory.  Finally, 
the  solid  curve  with  two  dots  represents  a  stiffened  shell  with  a  single 
equivalent  Poisson's  ratio  for  extension  ( =  0.  115)  used  in  both 
layers  (Ref.  11)  and  is  from  14  to  18  percent  lower  than  the  results 
of  the  present  theory.  The  lower  results  for  and  are  due  to 
neglect  of  coupling  between  bending  and  extension  of  the  two  shell 
layers.  Note  that  ai.  approaches  previous  to  the  present  theory  are 
conservative  for  this  jxample,  i  c.,  they  yield  lower  buckling 
pressures  than  can  actually  be  realized  by  the  stiffened  shell.  For 
other  problems,  the  previous  approaches  can  yield  unconservative 
results  (Ref.  11).  Thus,  the  importance  of  coupling  between  bending 
and  extension  should  not  be  overlooked. 
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SECTION  IV 


CONCLUDING  REMARKS 

An  exact  solution,  within  the  framework  of  classical  stability  theory, 
is  derived  for  the  buckling  of  a  circular  cylindrical  shell  with  multiple 
orthotropic  layers  and  eccentric  stiffeners  under  axial  compression, 
lateral  pressure,  or  any  combination  thereof.  The  simply  supported  edge 
boundary  conditions  are  6Nx  =  v  =  w  =  6Mx  =  0.  Thus,  the  present 
solution  can  be  regarded  as  a  lower  bound  on  results  for  practical  shells 
if  initial  imperfections,  prebuckling  deformations,  and  effects  of  discrete 
stiffener  spacing  are  ignored. 

A  numerical  example  is  given  to  illustrate  the  effect  of  coupling 
between  bending  and  extension  due  to  the  presence  of  different  layers  in 
the  shell  and  to  the  presence  of  eccentric  stiffeners.  Comparison  of  the 
present  theory  is  made  with  previous  approaches  such  ns  use  of  a  single 
equivalent  Poisson's  ratio  in  all  layers  of  a  layered  shell  and  orthotropic 
treatment  of  stiffened  shells.  The  buckling  predictions  of  the  previous 
approaches,  in  which  coupling  is  neglected,  ai  e  seen  to  be  erratic  in  that 
they  are  sometimes  conservative  and  sometimes  unconservative.  Thus, 
the  importance  of  coupling  between  bending  and  extension  should  not  be 
overlooked. 
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APPENDIX  A 


DESCRIPTION  OF  COMPUTER  PROGRAM 


A  computer  program  was  written  to  evaluate  the  closed-form 
stability  criterion,  Eq.  (18),  for  an  arbitrary  range  of  values  of  the 
buckling  mode  parameters  m  and  n  and  to  select  subsequently 
the  lowest  buckling  load  in  the  range.  Program  card  decks  are 
available  upon  request  to  the  Aerospace  Corporation,  San  Bernardino 
Operations,  Mathematics  and  Computation  Center.  Specific  charac¬ 
teristics  and  the  usage  of  the  program  are  described  in  the  following 
discussion. 

A.  1  GENERAL  CHARACTERISTICS 

The  basic  capability  of  the  program  is  represented  by  Eq.  (18) 
which  is  valid  for  the  stability  of  circular  cylindrical  shells  with 
multiple  orthotropic  layers  and  eccentric  stiffeners  under  axial 
compre ssion,  lateral  pressure,  or  hydrostatic  pressure.  The 
boundary  conditions  at  the  edges  are  6  N  =  v  =  w  =  6  M  =  0.  The 

/C  X 

orthotropic  material  properties  for  each  layer  of  thickness,  tk  , 

are  Exx  ’  Eyy  ’  \y  ’  Vyx  (reca11  that  because  of  the  reciprocal 
relations  only  three  are  independent)  and  Gky  .  It  should  be  noted 
that  the  principal  axes  of  orthotropy  must  coincide  with  the  shell 
coordinates.  The  geometrical  properties  for  the  stiffeners  are: 
area  (A),  moment  of  inertia  about  the  stiffener  centroid  (I),  eccen¬ 
tricity  (z),  torsional  constant  (J)  ,  and  spacing.  The  stiffeners  are 
isotropic,  hence,  E  and  v  are  the  only  material  properties  required. 

Because  mainly  algebraic  operations  are  performed  in  the  pro¬ 
gram,  the  execution  time  is  very  small  (less  than  1  second  per  case). 
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As  far  as  is  possible,  mnemonic  representations  are  used  throughout 
the  program. 

A.  2  ORTHOTROPIC  STIFFNESS  LAYER,  QSL 

Block,  Card,  and  Mikulas  included  an  orthotropic  stiffness 
layer  in  their  theory  (Ref.  3)  in  order  to  treat  corrugated  shells,  etc. 
In  the  present  program,  a  similar  layer  can  be  used  in  place  of  the 
first  layer  of  the  multilayered  shell  if  the  reference  surface  is  chosen 
to  be  the  middle  surface  of  the  orthotropic  stiffness  layer.  The 
orthotropic  stiffness  definitions  reduce  to  the  usual  definitions  for 
an  isotropic  shell,  1.  e. , 


Bx  =  By  =  B  =  Et/(  1  -  V2) 

B  =  [(1  -  v)/2]  B  =  Et/[2(1  +  v)] 
Dx  =  Dy  =  D  =  Et3/  [12(1  -  v2)] 

D '  =  [d  -  v)/2]  D  =  Et3/[24(1  +  v )] 

xy 


> 


(A- 


vxyB  vyxB  VxyD  WyxD  V  J 

The  orthotropic  stiffnesses  must  satisfy  the  reciprocal  relations 

VxyBBx  =  VyxBBy  and  VxyDDx  =  VyxDDy  *  K  iS  imPortant  to  note 
that  f  etc  are,  in  some  cases,  not  solely  material  properties, 

but  are  also  affected  by  the  geometry,  e.  g.  ,  corrugated  or  layered 

shells. 

The  orthotropic  stiffness  layer  was  used  to  describe  the  two¬ 
layered  eccentrically  stiffened  shell  in  Section  III,  Numerical 
Example,  in  order  to  obtain  the  curve  labeled  Orthotropic  Stiffness 
Approach  in  Figure  3.  Note  that  this  approach  neglects  coupling 


between  bending  and  extension  of  the  stiffeners  and  the  layered  shell 
and  also  neglects  coupling  between  bending  and  extension  of  the 
layers. 

Eccentric  stiffeners  can  be  added  to  the  orthotropic  stiffness 
layer  if  the  eccentricity  is  properly  accounted  for.  The  eccentricity, 
ZR  or  ZS  ,  is  ordinarily  input  as  the  distance  from  the  centroid  to 
the  base  of  the  stiffener.  Subsequently,  the  eccentricity  is  adjusted 
in  the  program  to  be  the  distance  from  the  centroid  of  the  stiffener  to 
the  arbitrary  reference  surface  of  the  layered  shell.  However,  when 
the  orthotropic  stiffness  layer  (0SL)  is  used,  the  reference  surface 
is  fixed  at  the  middle  surface  of  the  f£)SL  .  In  order  that  the  stiffener 
bend  about  the  middle  surface  of  the  layer  to  which  it  is  attached, 
it  Is  necessary  to  modify  the  input  eccentricity  such  that,  when 
one-half  the  (£)SL  thickness  is  added,  the  eccentricity  totals  one- 
half  the  thickness  of  the  layer  to  which  it  is  attached  plus  the  distance 
from  the  base  to  the  centroid  of  the  stiffener. 
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3  INPUT  PARAMETERS 


The  following  is  a  list  of  input  parameters  and  their  format  and 
definitions: 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

'  c 
c 
c 
c 
c 
c 


•CARO  1  FORMAT (BOH  I  -  PROBLEM  TITLE 
•CARO  2  FORMAT ( 1 10, 6F 10.01 

HL  -  NUMBER  OF  LAVERS  INCLUOING  ORTHOTROPIC  STIFFNESS  LAVER 
•RESTRICTED  TO  4  IN  0JMENS10N  LNI4)  AND  BV  FORMAT  NO.t.  THE  USUAL  THIN 
SHELL  LIMITATIONS  MUST  BE  TAKEN  INTO  CONSIDERATION  AS  NELL. 

OSL  -ORTHOTROFIC  STIFFNESS  LAVER 
IF  EQUAL  TO  0. ,N0  OSL 
IF  EQUAL  TO  1..0SL  REPLACES  LAVER  ONE 
LOAD  -  COCE  NAME  FOR  TVFE  OF  LOAO 
IF  EQUAL  TO  I.,  AXIAL  COMPRESSION 
IF  EQUAL  TO  2. <  LATERAL  PRESSURE 
IF  EQUAL  TO  3.,  HYDROSTATIC  PRESSURE 
NOfNF  -  INITIAL  ANO  FINAL  VALUES  OF  Mi  THE  NUMBER  OF  AXIAL  HALF-HAVES 
•HO  CANNOT  BE  ZERO  IN  THE  AXIAL  ANO  HYDROSTATIC  LOADING  CONDITIONS. 

•MO  SHOULD  BE  1  FOR  FINITE  LENGTH  SHELLS. 

•IF  NO  ABSOLUTE  MINIMUM  LOAO  IS  FOUND  OR  IF  THE  RELATIVE  MINIMA  ARE 
DECREASING  WHEN  N-MF,  A  MESSAGE  IS  PRINTED  STATING  THAT  TKE  RANGE 
ON  M  IS  INSUFFICIENT  TO  DETERMINE  AN  ABSOLUTE  MINIMUM. 

•THE  INTERVAL  CM0,4*MF>  IS  EXAMINED  INDEPENDENTLY  FOR  THE  AXISVNMETRIC 
BUCKLING  LOAO  HHICH  IS  THEN  PRINTED  ANO  ALSO  SAVED  FOR  COMPARISON 
WITH  THE  ASYMMETRIC  BUCKLING  LOAO. 

•THE  LONGER  THE  SHELL,  THE  hIGHER  MF  MUST  BE. 

NO.NF  -  INITIAL  ANO  FINAL  VALUES  OF  N,  THE  NUMBER  OF  CIRCUMFERENTIAL 
HAVES 

•THE  ENTIRE  INTERVAL  (NO.NF)  IS  EXAMINED  EVEN  IF  A  RELATIVE 
MINIMUM  IS  FOUNO  MITHIN  THE  INTERVAL. 

•NO  IS  NORMALLY  2  BECAUSE  A  SEARCH  FOR  THE  AXISVNMETRIC 
BUCKLING  LOAD  IS  AUTOMATICALLY  PROVIDED  IN  THE  AXIAL 
ANO  HYDROSTATIC  PRESSURE  LOADING  CONDITIONS. 

•NO  CANNOT  BE  ZERO  IN  THE  LATERAL  PRESSURE  LOADING  CONDITION. 
•MO  ANO  NO  CANNOT  BOTH  BE  ZERO  IN  THE  HYDROSTATIC  PRESSURE 
LOADING  CONDITION. 


.  C  • IF  NO  RELATIVE  MINIMUM  IS  FOUNO  OR  THE  LOAD  IS  AGAIN 

*  C  DECREASING  AFTER  ONE  MINIMUM  HAS  BEEN  FOUND  WHEN  N-NF, 

C  A  MESSAGE  IS  PRINTEO  STATING  THAT  TKE  INTEVAL  IS  INADEQUATE. 

C  *THE  THINNER  THE  SHELL,  THE  HIGHER  NF  MUST  BE. 

C 

C  *CARD$  3  THROUGH  NL»2  -  FORMAT UE  10.3 )  -  ORTHOTROPIC  LAVER  PROPERTIES 
C  LNll)  -  LAVER  NUMBER 

C  EXXI1I  -  MODULUS  OF  ELASTICITY  OF  THE  ITH  LAYER  IN  THE  X-OIRECTION 

C  EYYII)  -  MODULUS  OF  ELASTICITY  OF  THE  ITH  LAYER  IN  THE  Y-OIRECTION 

C  NUXYOI  -  POISSON'S  RATIO  FOR  CONTRACTION  IN  THE  Y-OIRECTION  DUE  TO 

C  TENSION  IN  THE  X-OIRECTION 

C  NUYXIII  -  POISSON'S  RATIO  FCR  CONTRACTION  IN  THE  X-OIRECTION  DUE  TO 

C  TENSION  IN  THE  Y-OIRECTION 

C  ‘NOTE  THAT  BY  THE  RECIPROCAL  RELATIONS  NUXV*EXX-NUVX*EYV. 

C  CXVIII  -  SHEAR  MODULUS  OF  ITH  LAVER  FOR  THE  XV-®LANE. 

C  TUI  -  THICKNESS  OF  THE  ITH  LAYER 

C  *IF  AN  ORTHOTROPIC  STIFFNESS  LAYER  IS  USED,  ALL  PROPERTIES  OF  THE 

C  FIRST  LAYER  ARE  ZERO. 


22 


ooooonnonooooooooooonno oononnnnooooo 


•CARO  0SL«<NL«3>  -  FORMATI5E10.3)  -  ORTHO  TROPIC  STIFFNESS  LAYER  FROFERTtES 
tx  -  EXTENSIONAL  STIFFNESS  IN  X-DIRECTION 
•Y  -  EXTENSIONAL  STIFFNESS  IN  V-DIRECTION 
•xv  —  <hC|i  tN  XY*FLAMf 

NUXYB-  EXTENSIONAL  POISSON'S  RATIO  FOR  CONTRACTION  IN  THE  V-OIRECTIQM 
OUE  TO  TENSION  IN  THE  X-OIRECTION. 

TOSL  -  MAXIMUM  THICKNESS  OF  OSL  (USED  AS  Tilt  IN  STIFFNESS  EQUATIONS 
FOR  LAVEREO  CYLINOER) 

•CARO  OSL*(NLMI  -  FORMAT  (4E10.3)  -  OSL  PROPERTIES,  CONTINUED 
OX  -  SENDING  STIFFNESS  IN  X-OIRECTION 
OV  -  SENDING  STIFFNESS  IN  T-OIRECTION 
OXV  -  TWISTING  STIFFNESS  OF  XY-FLANE 

NUXYD-  SENDING  PCISSON'S  RATIO  FOR  CURVATURE  IN  THE  V-OIRECTION 
OUE  TO  MOMENT  IN  THE  X-OIRECTION 

•CARO  NL*2*0SL*3  -  FORMATI6EI0.3)  -  RING  FROFERTIES 
ER  -  HOOULUS  OF  ELASTICITY 
AR  -  CROSS-SECTIONAL  AREA 

ZR  -  ECCENTRICITY  (NEASUREC  NEGATIVELY  INWARD  FROM  INNER  SURFACE  OF 
COMPOSITE  SHELL  TO  RING  CENTROID  IF  RINGS  ARE  INTERNAL  - 
POSITIVELY  OUTWARD  FROM  OUTER  SURFACE  IF  RINGS  ARE  EXTERNAL) 

IR  -  MOMENT  OF  INERTIA  OF  RING  ABOUT  ITS  OWN  CENTROID 
GRJR-  SHEAR  MOOULUS*TORS ION  CONSTANT  OF  CROSS  SECTION 
A  -  SPACING  OF  RINGS 

•CARO  NL*2*OSL*4  -  FORMAT (6E10. 3)  -  STRINGER  FROFERTIES 

EStASt ZS« IS.GSJStB  -  CORRESPCND  TO  ABOVE  RING  PROPERTIES 


•CARD  NL*2*OSl*5  -  FORMATI3E10.3I  -  BASIC  6E0METRY 
L  -  LENGTH  OF  CIRCULAR  CYLINDRICAL  SHELL 

R  -  RADIUS  TO  REFERENCE  SURFACE 

•MUST  BE  TO  MIOOLE  SURFACE  OF  OSL  IF  AN  OSL  IS  PRESENT 
DELTA-  DISTANCE  FROM  INNER  SURFACE  OF  LAYERED  CYLINDER  TO  REFERENCE 
SURFACE 

•MUST  BE  1/2*0 SL  THICKNESS  IF  AN  OSL  PRESENT. 

•SHOULO  GET  DIFFERENT  AXIAL  BUCKLING  LOAOS  WHEN  DELTA  VARIEO. 


A.  4  OUTPUT 


The  output  for  each  case  is  printed  on  one  page  if  the  sum  of 
the  number  of  layers,  LN,  and  the  number  of  axial  buckle  halfwaves, 
M  ,  does  not  exceed  25  and,  if,  in  addition,  there  is  no  more  than 
one  relative  minimum  buckling  load  per  value  of  M  .  If  these  con¬ 
ditions  are  not  met,  additional  pages  are  used  as  needed. 

First,  a  user- specified  case  identification  is  printed.  Next,  the 
input  quantities  are  printed  so  that  input  errors  can  be  identified. 

The  orthotropic  layer  properties  are  printed  and  are  followed  by  the 
orthotropic  stiffness  layer  (£)SL)  properties,  if  any.  Next,  the 
ring  and  stringer  properties  are  printed.  Finally,  the  basic 
geometry  quantities,  shell  length,  radius,  and  reference  surface 
location,  are  printed. 

After  execution  of  the  program,  the  buckling  load  for  axi- 
symmetric  deformation  (absolute  minimum  in  the  range  from  M  =  1 
to  M  =  4*MF)  is  printed  along  with  the  value  of  M  at  which  it 
occurs.  Subsequently,  the  asymmetric  buckling  loads  (relative 
minima  for  each  value  of  M  for  the  range  from  N  =  2  to  N  =  NF  ) 
are  printed.  The  final  result  is  the  absolute  rrrnimum  (axisymmetric 
or  asymmetric)  buckling  load  for  the  entire  range  of  M  and  N. 

A  typical  output  page  is  shown  ir.  Appendix  B. 


24 


APPENDIX  B 
EXAMPLE  PROBLEM 


The  example  chosen  here  is  the  configuration  discussed  in 
Section  III*  Numerical  Example,  in  the  main  body  of  the  report, 
i.  e. ,  a  ring- stiffened  circular  cylindrical  shell  with  two  isotropic 
layers  under  hydrostatic  pressure.  Pertinent  geometrical  and 
material  properties  are  given  in  Section  III.  Ring  spacing  for  this 
example  is  3  inches.  The  input  data  are  shown  in  Table  B-I. 
Figure  B- 1  illustrates  the  input  form,  and  the  computer  output  is 
shown  in  Figure  B-2. 
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Table  B-I 

INPUT  DATA  FOR  EXAMPLE  PROBLEM 


CASE  IDENTIFICATION: 
CONFIGURATION  OF 

FIGURE  3  -  ACTUAL  NUI 

Symbol 

Vasue 

Symbol 

i - 

Value 

NL 

0SL 

L0AD 

M0 

MF 

N0 

NF 

LN(1) 

EXX(l) 

EYY(l) 

NUXY(l) 

NUYX(l) 

GXY(l) 

T  (  1 ) 


44  x  10 
44  x  10* 
0 
0 

22  x  10* 
0.  04 


LN  (2) 
EXX(2 
EYY(2) 
NUXY  (2) 
NUYX(2) 
GX  Y  ( 2 ) 
T(2) 

ER 

AR 

ZR 

1R 

GRJR 

A 

L 

R 


DELTA 


2  x  10" 

2  x  106 

0.4 

0.4 

0.  7179  x  10* 
0.  3 

44  x  10^ 

0.  015 
-0.  125 
0.  7812  x  10' 
396 
3 

12 

6 

0.  02 


Example  Computer  Output 


APPENDIX  C 


FORTRAN  LISTING  OF  COMPUTER  PROGRAM 


c 

c 

c 

c 

c 

c 


c 


ELASTIC  BUCKLING  Of  SIMPLY  SUPPORTED,  ECCENTRICALLY  STIFFENED  CIRCULAR 
CYLINDRICAL  SHELLS  WITH  MULTIPLE  CRTHOTROPIC  LAYERS  UNDER  AXIAL  COMPRESSION, 
LATERAL  PRESSURE  OR  HYDROSTATIC  PRESSURE 


REAO  STATEMENT  FORMATS  -  BOLS 

1  Fl'RHAT  ( 8CH  BOLS 

l  BOLS 

2  FORMAT ( I 1C,7FL0.0!  BOLS 

3  l-ORMAT (8E1G.3)  6CLS 

WRITE  STATEMENT  FORMATS  -  BOLS 


4  FORMAT ( 90N  ELASTIC  BUCKLING  OF  SIMPLY  SUPPORTED,  ECCENTRICALLY  STIBOLS 
IFFENEO  CIRCULAR  CYLINDRICAL  SFELLS/57H  WITH  MULTIPLE  ORTHOTROPIC  LBOLS 


2AYERS  UNDER  AXIAL  COMPRESSION! 


BOLS 


5  FORMAT  1 9CH  ELASTIC  BUCKLING  OF  SIMPLY  SUPPORTED,  ECCENTRICALLY  STIBOLS 
IFFENEO  CIRCULAR  CYLINDRICAL  SFELLS/56H  WITH  MULTIPLE  ORTHOTROPIC  LBCLS 


2 AVERS  UNDER  LATERAL  PRESSURE! 


BOLS 


6  FORMAT (90H  ELASTIC  BUCKLING  OF  SIMPLY  SUPPORTED,  ECCENTRICALLY  STIBOLS 
1FF£‘  EO  CIRCULAR  CYLINDRICAL  SFELLS/60H  WITH  MULTIPLE  ORTHOTROPIC  LBOLS 


2AYERS  UNCER  HYDROSTATIC  PRESSURE!  BOLS 

7  FQRMATI/4H  MC-F4 .0 ,5X3HMF»F4.C , 5X3HNQ«F4.0»5X3HNF«F4.0 1  BCLS 

8  F0RMATC/15H  PROPERTIES  OF  ,Il,19H  ORTHOTROPIC  LAYERS/6H  LAYER, 7X3HB0LS 

1 EXX, X2X3HE YY , 12X4HNUXY , 1 1X4HNUYX , 1 IX3HGXY , 12XIHT i  BOLS 

9  F0RHATIF4. 0,4X513. 6, 5I2XE13. 611  BOLS 

10  FORMAT I/39H  ORTHOTROPIC  STIFFNESS  LAYER  PROPERTIES/SH  BX-E1 1. 4, 3X80LS 

13HBY«E 11 • 4, 3X4HB  XY«E 1 1 .4, 3X6HNUXYB*E 1 1 .4/5H  DX«E 11 .4, 3 X3H0V-E 1 1 .4B0LS 
2,3X4H0XY-Ell.4,3X6HNUXYD-EU.4,3X5HT0SL-EII.41  BOLS 

11  FORMAT (/ I6H  RING  PROPERT I ES, 32 X19HSTR I NGER  PROPERTIES/ 5H  ER-E11.4BCLS 

1,5X3H!R*EU.4,  15X3HES-E  1 1. 4,  5X3HI S«E  1 1 .4/5H  AR-E 1 1 . 4.3X5HGR  JR-E 1 1  BCLS 
2.4,1  5X3HA5*E 11 .4, 3X5HGSJS«E1 1. 4/5H  ZH^E 11 .4,6X2HA>E 11 .4 , 15X3HZS«EB0LS 
3U.4,6X2HB«E11.4!  BCLS 

12  FORMAT  ( /  15H  BASIC  C-E0HETRY/5H  L«E1 1.4,  3X2HR«E13.6,3X6H0ELTA*E12.BCLS 

15!  BOLS 

13  FORMAT I/22H  MINIMUM  NX  FCR  N»C  IS.E14.6.6H  AT  M»F4.0/9X1HM,7X21HREBCLS 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 
S6 

27 

28 
29 


29 


1LA7IVE  MINIM  »  NX.7X1MN)  SOLS 

14  FORMAT  I/21H  MINIMUM  P  FOM  M>0  1S.E14.4.4H  AT  N*F4.0//4XIMM,7X2OHRE*0LS 

1LATIVE  MINIMA  OF  P.tXIHN)  IOLS 

15  F0*NAT</4XlHMt7X20NRELATIVE  MINIMA  OF  P.IX1HN)  SOLS 

16  FORHAT(7XF4.0,6XE14.6«10XF4.0I  AOtS 

17  FOAM4TI/21H  At  SOLUTE  MINIMUM  NX-E14.4,5X2HH>F4.0«5X2HN>F4.0)  SOLS 

It  FORM  ATI /20H  ABSOLUTE  MINIMUM  PpE14.6,5X2HM»F4.0.SX2HI»F4.0)  tOLS 

C  EMOR  MESSAGE  FORMATS  tOLS 

14  FORMAT! 104H  THE  RELATIVE  MINIMA  ARE  STILL  OECREAS1NC,  SO  THE  RAN6EI0LS 
1  ON  M  IS  INSUFFICIENT  TO  OETERMINE  AN  AISOLUTE  MININUM/ltH  THE  LAStOLS 
2T  VALUE  IS.E14.6t6H  AT  M>F4.0)  tOLS 

20  FORMAT !62H  THE  LOAD  IS  DECREASING.  SO  THE  RANSE  ON  N  IS  INSUFFICIEtOLS 

1NT  TO  DETERMINE  ALL  MINIMA)  tOLS 

21  FORMAT  C/30H  EQUAL  OR  NEAR  EQUAL  0R01NATES/7XF4.0«tH  0R0NM1-E14.6,  BOLS 

110XF4.0/7XF4.0.BH  ORON-E 14.6, 10XF4.0// )  tOLS 

DIMENSION  Kill 4) .K12!4).K22<9) .K33I9) .UcLit) »EXX!4) »EVT!4I •NUXVItltOLS 
1.NUYXI9) .CXY!9).T!9)«LNi9)  tOLS 

REAL  IR» IS. M.HO. HF.MPL.N.NG.NF .NMtNUXY.MUYX.NUXVt.NUXVD.l.lOAO.  tOLS 
1K11.K12.K22.K33.LN.NM1  BOLS 

PI«3. 14159265  BOLS 

C  REAO  INPUT  DATA  BOLS 

100  REA0I5.1)  BOLS 

REAO(5.2)NL.OSL«LOA0.MO.MF.MO.kF  BOLS 

C  MRITE  TITLE  OF  DATA  AMO  PROBLEM  BOLS 

MRITE<6,1>  BOLS 

C  MRITE  TYPE  OF  LOAOING  ANO  RANGE  CM  N  ANO  N  BOLS 

IFJLOAO.EQ.l.l  NRITE<6.4>  BOLS 

lFtLGAO.EQ.2.!  WRITE(6f 5)  BOLS 

IFCLCAD.EQ.3.)  MR1TE(6.6)  BOLS 

MRITE! 6. 7)  HC.HF.NG.NF  BOLS 

C  REAO  ORTHOTROPIC  LAYER  PROPERTIES  BOLS 

00  110  1»1,NL  BOLS 

UC  RE  AD  ( 5 ,3 )  LMi).EXX<n,EYY(l>,  NUX  YU  I  .NUYX!  1 )  ,GXYC  I )  ,T(  !  )  BOLS 

IF(OSL.EQ.l..ANO.NL.EQ.l)  GO  TC  130  BCLS 

C  MRITE  QRTHOTRCP IC  LAYER  PROPERTIES  BOLS 

MRITE(6,6)  NL  BOLS 

00  120  I-l.NL  BOLS 

120  MR  IT  E!  6.  9)  LNI I  >  ,E  XX!  I  ),  EYY 1 1 )  ,NUX  V!  I J  .NUYX!  I )  ,GXY!  I ) ,  T  U)  BOLS 

C  TEST  FOR  PRESENCE  OF  CRTH0TR0P1C  STIFFNESS  LAYER  BCLS 

'r'n SI  T  ' .  '  Co  TO  13C  BOLS 

C  ZERO  OUT  PREVIOUS  ORTHOTROPIC  STIFFNESS  LAYER  PROPERTIES  BOLS 

BX-O.  BOLS 

BV»0.  BCLS 

8XY»0.  BOLS 

NUXYB=0.  BOLS 

TOSL’O.  BOLS 

0X»0.  BOLS 

0Y*O.  BCLS 

0XY*0.  BCLS 

NUXYD'O.  BOLS 

GO  TO  140  BCLS 

C  REAO  ORTHOTROPIC  STIFFNESS  LAYER  PROPERTIES  BCLS 

130  REAOI5.3)  BX,BY,BXY,NUXY8.T0SL  BOLS 

Till  *TOSL  BOLS 

RE  AO ( 5»  3 1  CX,OY,OXY,NUXYO  BOLS 

C  MRITE  ORTHOTROPIC  STIFFNESS  LAYER  PROPERTIES  BOLS 

MR  IT  El  6.  10 )  BX.BY.BXY.NUXYB.OX.OY.OXY.NUXYD.TOSL  BOLS 

C  REAO  ANO  MRITE  RING  AND  STRINGER  PROPERTIES  BOLS 

140  RE  AT  15*3)  £R  ,  AR,  ZR  ,  I  R.GR  JR  >  A  BOl  S 

RE AC! 5«3 1  EStASfZStISiGSJ5iB  80LS 

MRITE(6tlli  ERtIRiES,lS,AR,GRJR,AS(GSJS,ZR,A,ZS.B  BCLS 

C  REAO  ANC  MRITE  BASIC  GEOMETRY  BCLS 

REAO! 5(3)  L(Rt DELTA  BOLS 

MRITF.I6.12)  LtR. DELTA  *  BCLS 

C  CALCULATE  FUNCTIONS  OF  THE  ELASTICITY  CONSTANTS  BOLS 

DO  150  J  *  1 . NL  BOLS 

K1K  I  ) *£  XX !  I )  /  ( 1.— NUXY ( I  )*NUYX  ill)  BCLS 

K 1 21  I ) »NUX  VI I )  *K 1 1 ( 1 1  BOLS 


30 

31 

32 

33 

34 

35 
34 
37 
3t 

39 

40 

41 

42 

43 

45 

46 

47 
46 

49 

50 

51 

52 

53 

54 

55 

56 

57 
SB 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 
7! 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 


30 


R22I  D-EVYI 1 1/ 1  l.-NUXYI  1  )*NUY1<  1 1  > 

SOLS 

98 

190  x33<n*cxr<i> 

BOLS 

99 

CALCULATE  0EL*S  OF  THE  VARIOUS  LA TEAS 

BOLS 

100 

00  160  I«1»NL 

SOLS 

101 

IFII.EQ.I)  OEL(l)-m) 

BOLS 

102 

1FI1.NE.U  OH  II1-0EII  1-11*711 1 

BOLS 

103 

160  CONTINUE 

BOLS 

106 

ADJUST  ZR  AM  2$  TO  REFERENCE  SURFACE 

SOLS 

109 

IF1ZR.CT.0.I  ZR  -  ZR*(  0£L<  Nll-DELTA) 

BOLS 

106 

IF1ZS.GT.0.I  ZS  -  ZS*(OELC HL 1-0ELTA1 

BOLS 

107 

1FIZR.LT. 0.)  ZR  »  ZR-OELTA 

BOLS 

10B 

IF1ZS.LT.0.)  ZS  -  ZS-OELTA 

BOLS 

109 

CALCULATE  EXTENSIONAL,  COUPLING,  AMO  RENO INC  STIFFNESSES 

SOLS 

110 

ZERO  OUT  B»S,  C»S,  AND  D*S  PRIOR  TO  SUM NATION 

BOLS 

111 

SI 1-0. 

BOLS 

112 

B12-0. 

BOLS 

113 

822-0. 

BOLS 

116 

833-0. 

BOLS 

119 

Cl  1-0. 

BOLS 

116 

C12-0. 

BOLS 

117 

C22-0. 

BOLS 

118 

C33-0. 

BOLS 

119 

011-0. 

BOLS 

120 

012-0. 

BOLS 

121 

022-0. 

BOLS 

122 

033-0. 

BOLS 

123 

00  190  I  -  l.NL 

BOLS 

126 

IFII.NE.il  CO  TO  170 

BOLS 

125 

EXT-OELl 11 

BOLS 

126 

COUP- 1 1 . /2 . 1 • ( DEL ( 1 1 **2- 2. *DEL T A*OEL 1 1 t ) 

BOLS 

127 

8EN0-1 1./3. I*(0€L( U**3-3.*0EL1A*0EL ( ll**2*3.*DELTA*62*0ELClll  BOLS 

128 

CO  TO  180 

acts 

129 

170  EXT-OELl  Il-CEUI-l) 

BOLS 

130 

COUP-1  1./2.  )•!  (DELI I1P-2-0EL11-1  )**2)-2.*0ELTA*IDELI  I)-OEL(I-l))  )  BOLS 

131 

BEND-1  X./3.)*(IDELU»**3-OEL(I-1)*»3)-3.*OELTA*<OEH  Il**2- 

DELI I-1I80LS 

132 

1**2JA3.*0EITA**2*(0£L1  II-OEH  I-1J1J 

BOLS 

133 

180  11 1-B1  l*Kl  *1  1  l*EXT 

BOLS 

136 

B12-B12«K121II*EXT 

SOLS 

135 

822-622*K22(II*EXT 

BOLS 

136 

B33-833+X32.1  1)*EXT 

BOLS 

137 

Cll-CX1*K11<II*C0UP 

BOLS 

13B 

C12»C12*Ki;<I)*COUP 

BOLS 

139 

C22«C22«K22  I)*COUP 

BOLS 

160 

C33-C33AK33  ll-COUP 

BOLS 

161 

011-01 l*Kl  l  11-BEND 

BOLS 

162 

012-012+K12' II-8EN0 

BOLS 

163 

022-022+X221 I }*8EN0 

BOLS 

166 

190  033-D33*K3 11 1 1-BEND 

BOLS 

165 

INITIALIZE 

BOLS 

166 

ABSMH-.7E35 

BOLS 

167 

IF1LUA0.E0.2.I  GO  TO  300 

BOLS 

168 

CALCULATE  AXI S YMNETR1 C  BUCKLING  LCAOS  UNDER  AXIAL  OR 

HYDROS 

TIC  BOLS 

169 

LOADING  FOR  A  RANGE  OF  MO  TO  6-MF ,  AND  PRINT  MINIMUM 

LOAO 

SOLS 

150 

INITIALIZE 

BOLS 

151 

AX IM-6. *MF 

BOLS 

152 

M-MO 

BOLS 

153 

0R0MMI-.8E35 

BOLS 

156 

0RrMMZ-.9E35 

BOLS 

155 

200  MPL-M*PI/L 

BOLS 

156 

CALCULATE  A  VALUES 

BOLS 

157 

An«tBmBX*ES*AS/6)*NPL*»2 

BOLS 

158 

A 12-0. 

BOLS 

159 

Al3«(B12*NUXYB*BX)*NPL/R*(Cll-ES*AS*Z$/Bl*NPL*-3 

BOLS 

160 

A22*PB33-8XV1*MPL-«2 

8CLS 

161 

A23-0. 

BOLS 

162 

A3  3*(OU*OX*eS-lS/B*ES-AS-ZS— 2/61-MPL— 6*t2./R>*C12*MPL— 2*IK/R*80LS 

163 

1*2)*(822*BY*ER-AR/A) 

SOLS 

166 

PART* A33*(  (A12-A23-A13-A221/IA  11-A22-A12  — 2)  1-A13+1 < A12-AI 

V  A 1 1*  A26CLS 

165 

* 


5 

j 

X 


t 

l 

l 


31 


131/1  AU*A22-A12**2)I*A23  MLS  UA 

TEST  EM  TTEE  Of  LOAO  INC.  CALCULATE  SUCKLING  LOAO  (Ml  OK  PRESSURE),  KOU  1*T 

AMO  STORE  LOAO  IN  AOORESS  OROM  IOROINATE  AT  ARSCISSA  III  SOLS  1*< 

IFtlOAD.KI.l.l  OROW«PART/MPL**2  MLS  1*9 

IF (LOAO. EQ. 3. I  CROM-P ART/( .5*MPL**2I  SOLS  170 

TEST  FOR  AO SOLUTE  MINIMUM  Alt  SYMMETRIC  AUCXt  ,NG  LOAO  IOLS  171 

OROMM^  IS  TME  OROINATE  AT  AOSCISSA  N-l  IOLS  172 

0R0MM2  THE  OROINATE  AT  AOSCISSA  N-2  IOLS  173 

TEST  TO  S U  WHETHER  OROM  IS  INCREASING  OR  OECREASINS  OOLS  17* 

IFIOROM.GT.ORONHII  60  TO  210  OOLS  173 

OROM  OECREASINS  FROM  OR  EQU*«L  TO  OR0MM1  OOLS  17* 

If (M.EO.AIIMi  WRITE!*. 19)  OROM.M  OOLS  ITT 

60  TO  230  OOLS  170 

OROM  INCREASIN6  FROM  0R0MM1  OOLS  173 

210  If (0R0NM2.6T .OROMN1I  60  TO  220  OOLS  P» 

NO  RELATIVE  MINIMUM  FOUND  IOLS  101 

SO  TO  230  IOLS  102 

TEST  FOR  AOSOl’ITE  MINIMUM  SOLS  103 

220  IF(0R0MM1.6T. AOSMINI  60  TO  230  OOLS  10* 

NEW  ABSOLUTE  MINIMUM  FOUND  OOLS  109 

ABSMIN-OROMMl  OOLS  10* 

ASSM-M-l.  OOLS  117 

A8SN-0.  SOLS  100 

230  If (M.EO.AIIMI  60  TO  2*0  SOLS  189 

STEF  M  SOLS  190 

M*H*1«  SOLS  191 

0R0MN2-0R0HH1  SOLS  192 

OROHMIOROH  SOLS  193 

SO  to  200  SOLS  19* 

WRITE  /iXISVKMETRIC  BUCKLING  LOAO  SOLS  195 

2*0  !F ILn*0. EQ.l . I  HRITE(6.13i  ASSRINf ASSN  SOLS  196 

IF(L0A'/.EQ.3.)  WRITE!*.  1*1  ABSRIN.ABSN  SOLS  197 

CALCULATE-  ASYMMETRIC  BUCKLIN6  LOADS  FOR  A  SPECIFIED  RAN6E  OF  N  AND  N  SOLS  198 

INITIAL  It  *  SOLS  199 

300  N« NO  SOLS  200 

AM0NN1". 5E35  BOLS  201 

IFIL0A0.EQ.2.)  WRITE!*. 15)  SOLS  202 

8E6IN  H  LOOf  SOLS  203 

310  MPL-N*PI/L  SOLS  20* 

INITIALIZE  FOR  N  LOOf  SOLS  205 

N*NQ  BOLS  206 

0RDNN1-.8E35  SOLS  20T 

0R0NM2o.9E35  SOLS  208 

BE61N  N  LOOf  SOLS  209 

320  NR«N/R  B(Hj  210 

CALCULATE  A  VALUES  SOLS  211 

A1 1»( Bll*BX»ES*AS/B)*HPl**2+( B33+SXY )*NR**2  BOLS  212 

Al2>(6l2«NUXV0*ax*633+8XY)*MPl*NR  BOLS  213 

A13-!612*NUXY8*8XI*NPL/R*<CU*ES*AS*ZS/B)*MPL**3MC12*2.*C33I*NPL*B0LS  21* 
INR**2  BOLS  215 

A22-!B33*8XY1*MPL**2»(822*8Y*ER*AR/A)*NR**2  BOLS  216 

A23»(C12*2.*C33)*NPL**2*NR*!B22»BV»ER*AR/AI*NR/R*(C22*ER*ARO2R/Al*80LS  217 
INR**3  SOLS  218 

A33»(Oll«0X+ES9IS/8»ES*AS*ZS**2/BI»NPL*T**<2.*(O12*NUXYO*0X)**.»(OBOLS  219 
133»OXV)*GSjS/B»G.'UR/AI*MPL**2*HR**2M022*DY«'ER*IR/A4ER»AR*ZR**2/A)BOLS  220 
2*NM****! 2. /R >*C12*NPL**2  +  ( i./R )*IC22*ERf LR*ZR/AS*NR**2M l./R**2)*l BOLS  221 

3822*BV*ER*AR/A I  BOLS  222 

PART»A33*1  !AI2*A23-A13*A22)/IAU*A22-A12**2M*A13+!(A12*A13-A11*A2B01S  223 

13)/(All*A22-A12**2))*A23  BOLS  22* 

TEST  FOR  TYPE  OF  LOAOING,  CALCULATE  BUCKLING  LOAD  INX  OR  PRESSURE),  BOLS  225 

AND  STORE  LOAO  IN  AOORESS  ORON  i OROINATE  AT  ABSCISSA  H)  BOLS  22* 

IFILOAO.EQ. l.I0R0H*PA»T/MPl**2  BOLS  22T 

IFILOAO.EQ.2.JORON-PART  R »NR* *2)  801S  22B 

IF (LOAD. EQ . J . I ORDN*PART/ IR*( . 5*HPL**2*NR**2) I  BOLS  229 

BEGIN  TEST  FOR  RELATIVE  MINIMA  ANG  ABSOLUTE  MINIMUM  BOLS  230 

UR0NM1  IS  THE  OROINATE  AT  ABSCISSA  N-l  BOLS  ’31 

CR0NM2  IS  THE  OROINATE  AT  ABSCISSA  N-2  BOLS  232 

TEST  FOR  EQUAL  OR  HEAR  EQUAL  ORDINATES  BOLS  233 


32 


C*  1 
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IF I ABSI 2  ,* (OMON-OMBNMI)/ CORORMCMONNl  1 1  .CT.«  IE-3 1  CO  TO  390  COLS  234 

OFOINATES  U[  CLOSC  ENOUCN  TO  CAUSE  TMMJCLE  IN  THE  SEARCH  FOR  ROLS  235 

RELATIVE  MINIM,  $0  REST  IM»»«*Timi  i$  TO  MRITE  ORRINATES  ROLS  234 

NMl-N-1.  ROLS  23T 

HRITEI 4,21)  R,0R8NNt,NM,N,0«n».N  ROLS  23R 

CO  TO  3R0  ROLS  235 

TEST  TO  SEE  WHETHER  ORON  IS  INCREASING  OR  0ECREAS1NC  ROLS  240 

330  IF(ORON.CT.ORONNl)  CO  TO  340  ROLS  241 

ORON  OE CREASING  ROLS  242 

IFIN.EC.NF)  WRITEI4.20)  ROLS  243 

CO  TO  3R0  ROLS  244 

ORON  INCREASING  ROLS  245 

340  IFIORONH2.CT.OAONN1)  CO  TO  350  ROLS  244 

NO  RELATIVE  MINIMUM  ROLS  247 

CO  TO  3R0  ROLS  24R 

TEST  FOR  ARSOLUTE  MINIMUM  ROLS  249 

AM0H1  IS  THE  ARSOLUTE  MINIMUM  VALUE  OF  ORON  IN  THE  NF-1  LOOM  ROLS  250 

350  IFiN.EO.MF-l.. ANO.ORONH1.LT. ARCNM1)  ANONM1*ORDNM1  ROLS  251 

IF IN. E4. MF. ANO.MO. ME. MF.ANO.ORONM1.lt. AM 0NM1)  WRITEI4.19)  ROLS  252 

360  IFI0R0NN1.CT .ARSHIN)  CO  TO  370  ROLS  253 

NEW  ABSOLUTE  MINIMUM  FOUND  ROLS  254 

ABSMIN-0R0NM1  ROLS  255 

ABSF.-N  ROLS  254 

ABSN-N-1.  ROLS  257 

3TC  RELM1N-CRDNM1  ROLS  259 

RELN*N-1.  ROLS  259 

WRITE  REIATIVE  MINIMUM  WITH  CORRE  3F0N0I NC  M  A  HO  N  ROLS  240 

WRITEI6, 16SN,RELMIN,RELN  ROLS  261 

3S0  IFIN.EQ.NF)  CO  TO  390  ROLS  262 

STEF  N  ROtS  263 
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APPENDIX  D 


BONDLESS,  LAYERED  SHELLS 


The  objective  is  to  define  a  mathematical  model  for  a  circular 
cylirarical  shell  of  multiple  isotropic  layers  with  no  bond  between 
the  layers.  This  configuration  is  of  interest  as  a  lower  bound  to 
layered  shells  with  shear- deformable  bonds  between  the  layers.  The 
Kirchlioff-Love  hypothesis  is  employed  in  all  previous  sections,  but  is 
valid  only  if  the  bonds  between  layers  are  non- shear-deformable. 
Accordingly,  certain  nev-  definitions  must  be  established.  It  is  con¬ 
venient  to  work  within  the  framework  of  the  orthotropic  stixfness 
layer  feature  of  the  computer  program  (see  Section  A.  Z  of 
Appendix  A).  Certain  stiffnesses  and  so-called  Poisson's  ratios 
must  be  defined,  namely,  quantities  associated  with  extension 
(Bx,  3^,  B  ,  and  v  g)  and  those  associated  with  bending 

(D  ,  D  D  and  v  —  )  . 

x’  y’  xy’  xyD 

The  extensional  stiffness  of  a  layered  shell  is  not  affected  by  the 
presence  or  absence  of  a  bond  between  the  layers,  i.  e.  ,  it  remains 

N 

(D-" 

k--l 


Similarly,  the  resistance  to  in-plane  shear  is  unaffected,  so 
N 

xy  " 


B  =  ?  B.  (1  -  v,  )  /2 
xy  L-4  k  k' 


(D-2) 


k=  1 
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th 

where  is  the  bending  stiffness  of  the  k  layer  about  its  own  middle 
surface.  Note  that  there  are  no  terms  such  as  occur  in  the  transfer  axis 
theorem  for  moments  of  inertia,  i.  e. ,  no  (area)  times  (distance  squared) 
terms.  Consequently,  the  bending  stiffness  is  greatly  decreased  from 
the  perfect  bond  case. 

The  consistent  definition  for  the  twisting  stiffness  follows  from 
the  stipulation  that  each  layer  independently  resists  twisting.  Thus, 

N 

D*y=]C  Dk<I-l'k)/2  ,D-7> 

k=  1 


In  analogy  to  the  situation  for  extension,  it  is  stipulated  that  the 
layers  have  the  same  changes  in  curvature,  i.  e. , 


k  =  1,  N 


(D-8) 


Then  the  so-called  Poisson's  ratio  for  bending  is  obtained  by  use  of 
the  moment-change  in  curvature  relations  as 


N 


Again,  as  with  v  ,  is  a  geometrical  as  well  as  a  material 

property. 

The  above  approach  implies  that  the  layers  have  the  same  dis¬ 
placements  and  the  same  curvatures,  i.  e. ,  all  layers  take  the  same 
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shape.  This  implication  is  reasonable  as  long  as  the  layers  do  not 
separate. 

When  the  layers  are  in  contact,  the  membrane  circumferen¬ 
tial  strain  is  essentially  the  same  in  all  layers  if  the  su. n  of  the  layer 
thicknesses  divided  by  the  radius  of  the  shell  reference  surface  small, 


i.  e. ,  a  thin,  layered  shell.  Thus,  under  lateral  pressure,  which  is 

carried  as  membrane  circumferential  stress,  t  ,  in  the  present 

til 

buckling  theory,  a  in  the  k  layer  is  proportional  to  the  exten- 

th 

sional  stiffness  of  the  k  layer.  Accordingly,  the  lateral  pressure 


on  each  layer  is  given  by 


(D- 10) 


where  p  is  the  lateral  pressure  on  the  layered  shell.  Thus,  as  a 
crude  lower  bound  to  the  case  of  a  bondless,  layered  shell,  each 
layer  must  be  thick  enough  to  resist  buckling  under  the  pressure 
determined  by  Eq.  (D-  10).  In  addition,  the  layered  shell  with  stiff¬ 
nesses  given  by  Eqs.  (D- 1),  (D-2),  (D-5),  (D-6),  (D-7),  (D-9)  must 
be  thick  enough  to  resist  buckling  under  p  . 

Eccentrically  stiffened,  bondless,  layered  shells  can  be  treated 
by  appending  stiffeners  to  the  orthotropic  stiffness  layer  in  the  manner 
discussed  at  the  end  of  Section  A.  2  in  Appendix  A. 


APPENDIX  E 


TWO-LAYERED,  BONDLESS  SHELLS  WITH  CIRCUMFERENTIAL 
CRACKS  IN  THE  OUTER  LAYER 


The  objective  is  to  define  a  mathematical  model  for  a  circular 
cylindrical  shell  which  has  two  unbonded,  orthotropic  layers  and  cir¬ 
cumferential  cracks  in  the  outer  layer  (see  Figure  E-l).  The  princi¬ 
pal  axes  of  orthotropy  must  coincide  with  the  shell  coordinate  axes. 
The  orthctropic  stiffness  layer  feature  of  the  computer  program  (see 
Section  A.  2  of  Appendix  A)  is  used  in  the  calculations.  Accordingly, 
certain  stiffnesses  and  so-called  Poisson's  ratios  must  be  defined, 

namely,  quantities  associated  with  extension  (B  ,  B  .  B  ,  and 

x*  y’  xy* 

v„rR)  anc*  th°se  associated  with  bending  (D  ,  D  ,  D  ,  and  v  _). 
xyti  °  '  X*  y*  xy’  xyD7 

Because  of  the  circumferential  cracks  in  the  outer  layer,  the 

extensional  stiffness  in  the  axial  direction  is  merely  that  of  the 

inner  layer,  i.  e.  .  =  0.  However,  both  layers  are  effective  in 

resisting  circumferential  extension.  Thus, 

B  =  B  . 
x  xl 

(E-l) 


B  =  B  .  +  B  , 

y  yi  y2 

No  axial  strain  develops  in  the  outer  layer,  i.  e.,  t  =  0.  Thus, 
the  force-strain  relations  are 

Nx  Bxl  ^£xl  +  vxyBl  £yl^ 

Ny  "  Byl  (<yl  +  vyxB1  ‘*1>  +  By2  ‘  y2 


(E-2) 
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Moreover,  because  the  layers  do  not  separate  circumferentially. 


yl  y2  y 

Accordingly,  the  force -strain  relations  become 
N  =  B  ( €  +  v  e  ) 

x  x  x  xyB  y' 

N  =  B  (t  +  v  t  \ 

y  y  y  yxB  x' 

where  and  By  are  defined  in  Eq.  (E-l),  and 


(E-3) 


(E-4) 


vxyB  "  vxyBl 


VB  =  vyxBi  Byl  /(Byl  +  By2' 
Note  that  the  reciprocal  relations 


(E-5) 


VxyB  VyxB  ®y 


(E-6) 


are  satisfied  for  the  two-layered  shell  because  they  are  satisfied 
for  the  inner  layer,  i.  e.  , 


VxyBl  Bxl  =  VyxB  1  Byl 


(E-7) 


For  an  isotropic  inner  layer,  Eq.  (E-7)  is  an  identity. 

The  inner  layer  carries  all  the  in-plane  shear  because  the 
outer  layer  is  cracked.  Thus, 


B  =  B  , 
xy  xy  1 


(E-8) 


For  an  isotropic  inner  layer, 
Bxyl  =  EjV2*1  +V  i> 


( E  -  9 ) 
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Reasoning  parallel  to  the  above  leads  to  the  following  defini¬ 
tions  for  the  quantities  associated  with  bending. 


D 

x 

D 

y 

vxyD 

VyxD 


D 


xl 


D  .  +  D  , 
yl  y2 


xyDl 


VyxDl  Dyl/(Dyl  +  ^2* 


xyl 


(E- 10) 


(E  -11) 


(E-12) 


where,  for  an  isotropic  inner  layer. 


Dxyl  El'l  /24<1+V 


(E-13) 


In  the  definitions  in  Eqs.  (E-10)  to  (E-12),  it  is  implicit  that 


X  ,  ^  X  =  X 
yi  y2  y 


(E-14) 


in  analogy  to  Eq.  (E-3).  Both  Eqs.  (E-3)  and  (E-14)  are  a  result  of 
no  circumferential  separation  of  layers.  In  addition,  it  should  be 
noted  that  the  bending  stiffnesses  of  the  layers  in  Eq.  (E-16)  are 
about  the  middle  surface  of  the  respective  layers  because  of  the  lack 
of  bonding  between  layers. 

Eccentrically  stiffened,  bondless,  layered  shells  with  circum¬ 
ferential  cracks  can  be  treated  by  appending  stiffeners  to  the  ortho¬ 
tropic  stiffness  layer  in  the  manner  discussed  at  the  end  of  Section 
A.  2  of  Appendix  A. 
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1.  Delete  Eq.  (D-10)  and  the  discussion  in  the  surrounding  paragraph  on 
page  38  as  the  shell  buckling  analysis  is  unduly  conservative  if  the 
deleted  considerations  are  imposed.  That  is,  an  inner  layer  when 
constrained  by  an  outer  layer  would  be  expected  to  buckle  at  a  very 
much  higher  load  than  that  of  the  unrestrained  shell  implied  by  the 
deleted  considerations.  The  buckling  load  of  the  constrained  shell 
would  be  expected  to  be  higher  than  that  determined  by  the  model 
discussed  in  Appendix  D.  Thus,  the  model  in  Appendix  D  appears 

to  be  the  most  reasonable  model  which  could  be  devised. 

2.  Replace  Appendix  E  (pages  39  through  42)  with  the  attached  revised 
pages. 
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APPENDIX  E 


TWO -LAYERED,  BONDLESS  SHELLS  WITH  CIRCUMFERENTIAL 
CRACKS  IN  THE  OUTER  LAYER 


The  objective  is  to  define  a  mathematical  model  for  a  circular  cylindrical 

shell  which  has  two  unbonded,  orthotropic  layers  and  circumferential  cracks 

in  the  outer  layer  (see  Figure  E-l).  The  principal  axes  of  orthotropy  must 

coincide  with  the  shell  coordinate  axes.  The  orthotropic  stiffness  layer 

feature  of  the  computer  program  (see  Section  A.  2  of  Appendix  A)  is  used  in 

the  calculations.  Accordingly,  certain  stiffnesses  and  so-called  Poisson's 

ratios  must  be  defined,  namely,  quantities  associated  with  extension  (Bx>  B^, 

B  ,  and  v  „)  and  those  associated  with  bending  (D  ,  D  ,  D  ,  and  v  ~). 
xy  xyB'  &  '  x  y  xy  xy D' 

Because  of  the  circumferential  cracks  in  the  outer  layer  and  the  lack 
of  bonds  between  layers,  the  axial  force  in  the  outer  layer  is  zero,  i.e.  , 


Nx2  =  Bx2  (€x2  +  VxyB2  ey2)  =  0 


(E-l) 


The  remaining  segments  of  the  outer  layer  are  analogous  to  plane  stress 
ring  elements,  the  axial  stiffness  of  which  is  finite.  Accordingly,  from 
Eq,  (E-l), 


x2 


VxyB  2  6  y2 


(E-2) 


The  force -strain  relations  can  then  be  written  as 


N  =  B  .  (e  .  +  v  „.«.) 
x  xl  xl  xyBl  y  1 


N  =  B  .  (c  .  +  V  Tal€.)  +  B-(€,1-V  „,  «  -) 

y  yl  yl  yxB  1  xl  y2  y2  yxB2  x2 

Moreover,  because  the  layers  do  not  separate  circumferentially, 


(S-3) 


£  ,  =  €  ,  =  C 

yl  y2  y 


(E-4) 
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whereupon,  with  Eq.  (E-2),  the  force  strain  relations  become 

N  =  B  (e  +  v  x,  e  ) 
x  x  x  xyB  y 

N  =  B  («  +  v  „  €  ) 
v  y  y  yxB  x 


(E-5) 


where 


B  =  B  . 
x  xl 


By  ~  Byl  +  By2  {1  '  VyxB2  VxyB2} 


(E-6) 


€x  =  cxl 


VxyB  “  VxyBl 

v  _  =  v  „ .  B  /B 
yxB  yxB  1  yl  y 

Note  that  the  reciprocal  relations 

v  „  B  =  v  B 

xyB  x  yxB  y 


(E-7) 


(E-8) 


are  satisfied  for  the  two-layered  shell  because  they  are  satisfied  for  the 
inner  layer,  i.  e.  , 


v  _ .  B  =  v  „ .  B  , 
xyBl  xl  yxBl  yl 

For  an  isotropic  inner  layer,  Eq.  (E-9)  ie  an  identity. 


(E-9) 


The  inner  layer  carries  all  the  in-plane  shear  because  the  outer  layer 
is  cracked.  Thus, 


B  =  B  . 
xy  xyl 

For  an  isotropic  inner  layer, 
Bxyl  =ElV2<l  +  vl> 


(E  - 10) 


(E  -11) 


41 


* 

4 


Q 

LU 


LU 


Reasoning  parallel  to  th^  above  leads  to  the  following  definitions  '65 
the  quantities  associated  with  bending:  ^ 


D  =  D  . 
x  xl 


Dy  "  Dyl  +  Dy2  (1  '  VyxD2  VxyD2) 


(E-12) 


D  =  D  , 
xy  xyl 


and 


vxyD  VxyD  1 

v  r.  =  v  ,  D  ./D 
yxD  yxDl  yl  y 

where,  for  an  isotropic  layer, 

Dxyl  =  E1  24“  +  »l> 

In  the  definitions  in  Eqs.  (E-12)  and  (E-13),  it  is  implicit  that 
*yl  =  *y2  =  *y 


(E-13) 


(E- 14) 


(E- 15) 


and 

=  Xx  (E-16) 

in  analogy  to  Eqs.  (E-4)  and  (E-6).  Both  Eqs,  (E-4)  and  (E-15) 
are  a  result  of  no  circumferential  separation  of  layers,  in  addition,  it 
should  be  noted  that  the  bending  stiffnesses  of  the  layers  in  Eq.  (E-12) 
are  about  the  middle  surface  of  the  respective  layers  because  of  the  lack 
of  bonding  between  layers. 

Eccentrically  stiffened,  bondless,  layered  shells  with  circumferential 
cracks  can  be  treated  by  appending  stiffeners  to  the  orthotropic  stiffness 
layer  in  the  manner  discussed  at  the  end  of  Section  A.  2  of  Appendix  A. 
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